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Abstract. This paper presents new near-infrared obser- 
vations of the planetary nebula NGC 2346. The data in- 
clude a broad K-band image, an image in the H2 vibra- 
tional!}^ excited 1-0S(1) line and K band slit spectra at 
three positions in the nebula. In the H2 1-0S(1) line, the 
nebula is characterized by a central, bright torus, sur- 
rounded by weaker emission with a typical butterfly shape, 
as seen in Ha and CO lines. The K band spectra show 11 
H2 lines with excitation energies from 6150 to 12552 K. 
The H2 data have been compared to the predictions of 
models which follow the evolution with time of the H2 
emission in PNe of different core mass and shell proper- 
ties (Natta & Hollenbach 1998). These models compute 
the emission originating in the photodissociation region 
(PDR) created at the inner edge of the neutral shell by 
the UV radiation of the central core, as well as the emis- 
sion in the shock associated to the expansion of the shell 
inside the precursor red-giant wind. In NGC 2346, a PDR 
origin of the H2 emission in a low-density molecular shell 
(n < 10 4 cm~ 3 ) is indicated. At these low densities, time- 
dependent H2 chemistry and X-ray heating of the neutral 
gas enhance the predicted PDR H2 line intensity by large 
factors. 

Key words: Planetary Nebulae; NGC 2346; IR spec- 
troscopy; H2 lines 



1. Introduction 

NGC 2346 is a much studied bipolar planetary nebula at a 
distance D ~ 800 pc (Acker et al. 1992)0. At its center lies 

Send offprint requests to: natta@arcetri.astro.it 
* Based on observations obtained at the TIRGO telescope 
1 After this work was completed, a new determination of 

the distance (_D=690 pc) was quoted by Terzian (1997). It is 

close enough to the value we adopt in this paper that none of 

the conclusions needed to be changed. 



a binary system, formed by a main-sequence star of spec- 
tral type A5V, with mass ~ 1.8 M©, temperature ~ 8000 
K and luminosity ~ 14 L© (Mendez and Niemela 1981; 
Walsh 1983) and hot star, not detected in the visual, with 
T*~ 130000 K (Mendez 1978), which excites the nebula. 
Its luminosity is very uncertain, as we will discuss in §4. 
Estimates in the literature give L*— 17 — 90 L Q (Mendez 
1978; Calvet and Peimbert 1983). 

In the optical, the nebula has a butterfly shape (Balick 
1987; Walsh et al. 1991), with well developed bipolar lobes 
and a bright torus which surrounds the central star. The 
temperature and density of the ionized gas in the torus 
have been estimated to be T - 12000 K and n e < 700 
cm~ 3 , respectively (Liu et al. 1995; McKenna & Keenan 
1996). The nebula contains a large amount of material in 
the form of molecular g revealed by the CO obser- 
vations of Knapp (1986), Huggins & Healy (1986), Healy 
& Huggins (1988). Bachiller et al. (1989) have mapped 
the entire nebula in the two CO lines J=l-0 and J=2-l; 
the morphology of the molecular gas follows very closely 
that of the ionized gas, showing a clumpy, inhomogeneous 
torus, tilted with respect to the line of sight by about 56°, 
which is expanding outward. Scaled to our adopted dis- 
tance D = 800 pc, the torus has a radius of —0.05 pc, and 
mass ~ 0.26 M Q , much larger than the mass of ionized 
gas (-0.01 Mq, Walsh 1983). The radial velocity of the 
most intense CO condensations is of the order of 15-35 
km s^ 1 , which results in a dynamical age of about 2500 
yr (but see also Walsh et al. 1991). 

NGC 2346 is a Type I nebula, originated by a massive 
progenitor (Calvet and Peimbert 1983). As many PNe of 
similar type, NGC 2346 is detected in the vibrationally 
excited lines of H2 (Webster et al. 1988). Zuckerman &; 
Gatley (1988) have mapped the nebula in the 1-0S(1) line 
at 2.12 fim using a single-beam 12 arcsec spectrometer 
with resolution — 200. The morphology of the nebula in 
this line is again very similar to the morphology observed 
in the optical lines and in CO. This result was confirmed 
more recently by the images obtained in the same line 
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with much better spatial resolution (about 1-2 arcsec) by 
Kastner et al. (1994) and Latter et al. (1995). 

The excitation mechanism of the vibrationally excited 
H2 lines in this, as in other PNe, is still uncertain. Zucker- 
man & Gatley (1988) discuss the possibility that they form 
in a shock driven by the fast wind emitted by the central 
star. Kastner et al. (1994) surveyed a sample of bipolar 
planetary nebulae (including NGC 2346); they conclude 
that the H2 emission very likely originates in thermally 
excited (possibly shocked) molecular gas. Recently, Natta 
& Hollenbach (1998; hereafter NH98) have computed theo- 
retical models of the evolution of PN shells and predicted, 
among others, the intensity of the most commonly ob- 
served H2 vibrationally excited lines (namely, the 1-0S(1) 
at 2.12 /xm and the 2-lS(l) at 2.25/im). They consider the 
emission of the photodissociation region (PDR) formed by 
the UV photons emitted by the central star impinging on 
the shell, including in the calculations time-dependent H2 
chemistry and the effects of the soft X-ray radiation emit- 
ted by the central star, which are important in sources 
like NGC 2346 where T*>10 5 K. NH98 compute also the 
emission of the shocked gas at the interface between the 
shell and the wind ejected by the central star in its previ- 
ous red giant phase. They point out that both mechanisms 
(PDR and shocks) can produce lines of similar intensity, 
with reasonable values of the model parameters. 

The PN properties that determine the intensity of the 
H2 lines are very different in the two cases. As discussed 
in NH98, if the emission is produced in the warm, neutral 
PDR gas, the line intensity depends mostly on the stel- 
lar radiation field which reaches the shell and, to a lower 
degree, on the density of the neutral gas itself. If the emis- 
sion is produced in the shocked gas, then the line intensity 
does not depend directly on the properties of the central 
star or of the PN shell, but only on the shock velocity 
and on the rate of mass-loss of the precursor red-giant. It 
is therefore clear that, before attributing any diagnostic 
capability to the H2 lines, we need to understand which 
of the possible excitation mechanisms dominate the PN 
emission. 

This paper is a first attempt to understand the H2 
emission of a well-studied PN in a quantitative way, i.e., 
by comparing the observations to detailed models of PDR 
and shock emission, such as those discussed in NH98. To 
this purpose, we have collected new near-IR broad and 
narrow-band images of NGC 2346 as well as K band spec- 
tra with resolution ~ 1000. These observations are de- 
scribed in §2. The results are described in §3 and com- 
pared to the predictions of PDR and shock models in §4. 
A discussion of the results follows in §5; §6 summarizes 
the main conclusions of the paper. 



2. Observations 



2.1. ARNICA Observations 

NGC 2346 was observed during two observing runs in Jan- 
uary 1996 using ARNICA (ARcetri Near Infrared CAm- 
era) mounted on the 1.5m TIRGOQ telescope. ARNICA 
is equipped with a 256x256 NICMOS3 array, the pixel 
size with the optics used at TIRGO is 0.96"; for a com- 
plete description of instrument performances, see Lisi et 
al. Ql996j ) and Hunt et al. ( |1996j ). Ima ges were obtained 
in the K broad-band filter (centered at 2.2 /xm) and in 
a narrow-band filter cente red o n the 2.12 H2 1-0S(1) line 
(A A/ A ~ 1%, Vanzi et al. |l998| ). The seeing was approx- 
imately 2-3"and the observed field was ~ 2' x 2', covering 
all the nebula. 

Data reduction was carried out using the IRAF0 and 
ARNICA (Hunt et al. 1994) software packages. Photomet- 
ric calibration in the K band was performed by observing 
the photometric standards of the FS14 group from the 
list of Hunt et al. (1997). The quality of the night was 



rather poor, and the calibration accuracy is estimated to 
be ~ 15%. 

The image in the H2 1-0S(1) line has been calibrated 
using the 5 brightest (unsaturated) stars in the ARNICA 
images, under the assumption that for each star the flux 
density measured in the line filter was equal to the flux 
density measured in the K band. Integrated line fluxes 
on the nebula were then obtained multiplying the flux 
density by the bandwith of the narrowband filter (Vanzi 
et al. 1998). The accuracy is ~15%. 



2.2. LONGSP Observations 

K (2.2 [im) band spectra of NGC 2346 were obtained us- 
ing the LonGSp (Longslit Gornergrat Spectrometer) spec- 
trometer mounted at the Cassegrain focus on the TIRGO 
telescope. The spectrometer is equipped with cooled re- 
flective optics and grating in Littrow configuration. The 
detector is a 256x256 engineering grade NICMOS3 ar- 
ray (for detector performances see Vanzi et al. 1995). The 
pixel sizes are 11.5 A (first order) and l'/73 in the disper- 
sion and slit directions, respectively. LONGSP operates 
in the range 0.9-2.5 /im achieving a spectral resolution at 
first order of A/AA ~950 in K. For a more comprehensive 



description of the instrument, refer to Vanzi et al. (1997). 

Observations were conducted in March 1998 under non- 
photometric conditions. The slit had dimensions 3'.'5x70" 
and was oriented N-S. The seeing during the observations 
was in the range 2"-4". NGC 2346 was observed at three 



2 The TIRGO telescope is operated by the C.A.I. S.M.I.- 
C.N.R Firenze, Italy 

3 IRAF is made available to the astronomical community 
by the National Optical Astronomy Observatories, which are 
operated by AURA, Inc., under contract with the U.S. National 
Science Foundation 
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slit positions labeled as E, S, W and shown in Fig. [l] su- 
perimposed on the image in the H2 1-0S(1) line. Position 
E and W are centered on the peaks of the line emission 
located east and west of the central star, respectively (see 
Zuckerman and Gatley 1988). Position S is centered on 
the star. At each grating position we performed 5 ABBA 
cycles (A=on source, B=on sky) with an on-chip integra- 
tion time of 60 sec, for a total of 10 min integration on 
source. 

Data reduction was performed with the ESO package 
MIDAS, within the context IRSPEC, modified to take into 
account LonGSp instrumental characteristics. The frames 
were corrected for bad pixels, flat-fielded, sky subtracted 
and wavelength calibrated using the OH sky lines present 
in all the frames (Oliva & Origlia 1992). After direct sub- 
traction, sky removal was optimized by minimizing the 
standard deviation in selected areas where the OH sky 
lines were poorly subtracted but no object emission was 
present. The wavelength calibration was performed to bet- 
ter than 1/5 of a pixel (~2A) . The spectra were then cor- 
rected for telluric absorption by dividing the spectra by 
the spectrum of the A star BS 2714 after removing its 
photospheric features (mainly Br7). For more details on 
LonGSp data reduction, see Vanzi et al. 
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Flux calibration of the spectra was achieved by rescal- 
ing the observed flux distribution along the slit to match 
that obtained from the ARNICA image in the 1-0S(1) line 
at the positions of the slits. 




20 

Aa (arcsec) 

Fig. 1. Image in the H2 1-0S(1) line. The contours range 
from 1 to 14 in steps of 1.3 x 10~ 5 erg cm— 2 s _1 sr . 
The three positions of the slit used to obtain the K-band 
spectra are shown by the black boxes and are marked as 
E, S and W, respectively. The continuum (not subtracted 
from the image) is confined to a region of ~ 5" radius 
centered on the star. 



3. Results 

The image in the H 2 1-0S(1) line is shown in Fig. [l]. As 
verified with the spectra, the continuum emission is ev- 
erywhere negligible but at the position of the central star, 
i.e. in a 5" radius region centered on the star. 

The H2 1-0S(1) image of Fig. [I] shows the well-known 
NGC 2346 morphology, with a bright central region of 
size ~ 50"x20" and two very extended lobes of weaker 
emission (Kastner et al. 1994). The central region has two 
peaks of emission, to the east and west of the star and 
matches well the bright torus, tilted with respect to the 
line of sight, seen in optical tracers and in CO (Walsh 
1983; Bachiller et al. 1989). The H 2 1-0S(1) intensity is 
~ 1.3 x 10 -4 erg cm~ 2 s _1 sr — x on both peaks. The total 
luminosity of the nebula in this line is about 0.06 L/0(for 
D = 800 pc), of which about 48% is contributed by the 
torus. The average line intensity over the torus (defined as 
the central region of size 20" x 50") is ~ 6 x 10 -5 erg cm -2 
s _1 sr . These numbers are very similar (within 10-20%) 
to those derived by Zuckerman and Gatley (1988). 

The K band spectra in the three positions E, S and W 
are shown in Fig. |^. The spectra have been averaged over 
a region of ~20" (12 pix) along the slit centered on the 
torus midplane. The line intensities are given in Table 1 , 
which gives in Column 1 the line identification, in Column 



2 the wavelength of the line, in Column 3 the intensity 
in the W position, in Column 4 that on the E position, 
in Column 5 that in the S position of the slit. The lines 
are normalized to the 1-0S(1) line set equal to 100; the 
intensity of the 1-0S(1) line is given in the Table's note. 
Typical uncertainties on the line ratios are ~10% for ratios 
>50, and —30% for the others. Lines whose intensity is 
particularly uncertain are marked with a semicolon. 

The variations in intensity along the slit positions E 
and W and S of the H 2 1-0S(1) and 2-lS(l) lines and of 
Br7 is shown in Fig. [| The intensity profiles show that the 
H2 emission in the 1-0S(1) line peaks in the torus (i.e., for 
^1 — 10") an d has extended emission with intensity that 
declines more sharply toward the north than toward the 
south, possibly due to the tilt in the plane of the sky. 
Some of the inhomogeneities seen in the two-dimensional 
image can also be seen in these profiles. The ratio of the 
2-lS(l) over the 1-0S(1) intensity varies between about 
0.1 in the W and E peaks of emission to about 0.04 in the 
condensation detected in the W position —15" south of the 
midplane. Another interesting result which emerges from 
the spectra and the intensity profiles is that the emission 
in Br7 extends over the whole central torus. However, Br7 
is quite weak, -10% of the 1-0S(1). 
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Fig. 2. K-band spectra in the W, E and S positions. The 
lines have been averaged over a slit section of ±10" cen- 
tered on the peak of the emission (the torus mid-plane). 
The small gap at ~ 2-2. 04 /im corresponds to a region of 
very poor atmospheric transmission. 

4. Comparison of the H 2 spectrum with models 

The H 2 emission in the T0S(1) line of the central torus of 
NGC 2346 is compared in this section to the predictions 
of the NH98 theoretical models. 

The NH98 models compute the emission expected from 
tori (or shells or clumps) exposed to the radiation field of 
the hot central star and expanding inside the remnant 
of the wind ejected by the central star in its previous 
red giant phase. As the central star evolves initially at 
constant luminosity toward higher effective temperatures 
and then along the white dwarf cooling track, the torus 
expands radially with roughly constant velocity and de- 
creasing density. The models consider the effects of UV 
and soft X-ray radiation on the neutral gas and follow the 
time dependent chemistry for H2, solving for the chemical 
and temperature structure and the emergent spectrum (in 
the following, the PDR spectrum) of the evolving torus. 



Table 1. K-Band Spectrum of the Central Torus 
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Note: 1-0S(1) intensity 100 corresponds to 7 ± 2 x 10 5 erg 
cm _2 s _1 sr _1 in the W position, 9 ± 3 x 10~ 5 erg cm _2 s _1 
sr _1 in the E position and 3 ± 1 x 10~ 5 erg cm _2 s _1 sr _1 
in the S position. The spectra have been averaged over a slit 
portion ±10" centred on the torus midplane. 



The torus expands inside the material ejected by the 
star in its previous phase as red giant. Since its velocity 
is larger than that of the wind itself, it gives origin to 
a shock which heats and compresses the gas, which then 
emits intense lines of vibrationally excited H2. NH98 con- 
sider the emission of J shocks, since the importance of 
magnetic field in PNe is not clear and, in this range of 
shock velocities, the H 2 line intensity is maximum in J 
shocks. 

4.1. PDR Models 

The NH98 PDR models require the specification of a num- 
ber of parameters. We fix the mass of the central core, 
which determines the time scale of the evolution of the 
stellar radiation field, to be about 0.7M Q (Calvet and Pc- 
imbert 1983; Bachiller et al. 1989). For a core mass M*=0.7 
M Q , NH98 show the results of only one model, which 
has a neutral gas density at the time t=2500 yr (roughly 
the age of NGC 2346) n ~ 7 x 10 5 cnr 3 (note that in 
NH98 the density varies with time as n = n (t/t )~ 2 )- 
We used the NH98 code to compute a number of addi- 
tional models, varying the t — 2500 yr density no from 
7x 10 3 cm~ 3 (Model 1) to 2.1 x 10 4 cnT 3 (Model 2), 7x 10 4 
cm~ 3 (Model 3) and 7x 10 5 cm" 3 (Model 4). The model pa- 
rameters are summarized in Table 2. In all cases, He/H=0.13, 
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Fig. 3. Top panel: variation with declination offset of the 
intensity of the H 2 1-0S(1) and 2-lS(l) lines and of Br7 
(dotted curve) for the slit position E. Middle panel: same 
for slit position W. Bottom panel: same for slit position 
S. All intensity profiles have been smoothed over three 
pixels. The offset is measured with respect to the central 
star. 



C/H=5.3 x 1CT 4 , 0/H=4.6 x 1CT 4 . The results are shown 
in Fig. |J, which plots the intensity of the H2 1-0S(1) line 
as a function of time. 

The square in Fig. [| shows the observed intensity of 
the 1-0S(1) line averaged over the torus. We have esti- 
mated the age of NGC 2346 by taking the separation of 
the two 1-0S(1) peaks (about 30 ") at D = 800 pc and an 
expansion velocity of ~25 km s _1 (as in the NH98 models), 
consistent with the observed range (15-35 km s _1 ) of CO 
expansion velocity (Bachiller et al. 1989). The result is a 
dynamical age of ^2500 yr. The uncertainty on the age is 
certainly large, but difficult to estimate. We plot in Fig. [| 
an error bar corresponding to an uncertainty of ±500 yr. 
The observed intensity of the 1-0S(1) line is reproduced 
quite well by the low-density models, especially if we take 
into consideration the large uncertainties that affect the 
models as well as the age estimate of the nebula itself. 

One problem that arises immediately with these mod- 
els has to do with the estimated luminosity of the central 
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Fig. 4. Model predictions for the H2 T0S(1) line intensity 
as function of the age of the nebula. Model parameters are 
given in Table 2. All the models have a neutral gas den- 
sity decreasing as t~ 2 as the torus expands with constant 
velocity v — 25km s _1 . The square shows the value of the 
observed H2 1-0S(1) intensity averaged over the central 
torus. 



star. At t ~ 2500 yr, a core of 0.7 Mq has already reached 
the white dwarf cooling track; it has a luminosity of about 
250 L Q , an effective temperature T*^ 1.5 x 10 5 K and 
a number of ionizing photons 3 x 10 45 photons s _1 
(Blocker 1995). While the effective temperature is roughly 
in agreement with the Zanstra Hell temperature (Mendcz 
1978), the luminosity is significantly larger than the val- 
ues 17-90 L© quoted in the literature. However, we are 
somewhat suspicious of those very low values. The lumi- 
nosity derived from the Hell A4685A intensity by Mendez 
(1978) is a lower limit L*>43 L Q (for L>=800 pc) and is 
very sensitive to the extinction. The number of ionizing 
photons we derive from the observed radio flux at 6 cm 
(86 mJy; Milne and Aller 1975) and from the total Ha 
flux (Walsh 1983) is at least 2 - 4 x 10 45 photons s" 1 
(assuming no escape of ionizing photons and an average 
optical depth in Ha ~ 0.7), consistent with L*^250 L , 
but not with lower values of L*(see Fig. 3 of NH98). If 
L^=250 Lq, the non-detection of the white dwarf star in 
the visual is not surprising; assuming, for simplicity, that 
the A star and the white dwarf spectrum can be repre- 
sented by black-bodies at 8500 K and 1.5 x 10 5 K, having 
luminosities of 15 and 250 Lq, respectively, we find that 
the white dwarf is a factor 54 weaker than the A-type star 
at 5500 A, a factor of 10 at 3000A and that the two stars 
become comparable only at ^2000A. This last is consis- 
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Table 2. Model Parameters 
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(a): Mod 4 assumes a luminosity of the central core equal to 
L*(M*=0.7)/5 at all times. 



tent with the UV excess (with respect to the flux expected 
for the A star) measured by the ultraviolet satellite ANS 
and reported by Mendez (1978). All together, we suspect 
that the white dwarf luminosity is roughly of the order 
of 250 L/0. In any case, we have also computed a model 
where we have artificially reduced the stellar luminosity by 
a factor 5 at all times; the density of this model (Model 
5), shown as a dashed curve in Fig. [|, is no — 7 x 10 3 
cm -3 . The predicted line luminosity scales approximately 
with the luminosity of the central core. A value L*=50 
L Q (although not predicted by any evolutionary track) is 
still roughly consistent with the observed line intensity, 
especially if we consider the uncertainty on the PN age 
estimate. However, this model predicts an intensity lower 
than observed (by a factor 3-10) for all the H2 lines we 
measured. 

The best fit to the H2 1-0S(1) observations is provided 
by models with low density (no ~ 1 - 3 x 10 4 cm -3 ), 
in good agreement with the low electron density (<10 3 
cm -3 ) derived for the ionized part of the nebula (Liu et 
al. 1995; McKenna & Keenan 1996). Assuming pressure 
equilibrium between the ionized and the neutral gas, and 
a PDR temperature ~500 K, we expect a neutral density 
about 40 times the electron density. The low density we 
require is in rough agreement with the Bachiller et al. 
(1989) estimate that the neutral density is fewx 10 3 cm -3 . 

In these low-density models, at t ~ 2500 yr, the 1- 
0S(1) emission is mostly due to collisionally excited H2, 
kept warm by the heating of the soft X-rays emitted by the 
central core. As discussed in NH98, X-rays determine the 
chemical and physical evolution of the neutral gas around 
high-mass PN cores, after a short initial phase (about 1000 
yr for a 0.7 M core) where UV photons dominate. If the 
X-rays effects are neglected, PDR models predict a much 
lower intensity of the H2 molecular lines. 

Time-dependent effects in the H 2 chemistry are impor- 
tant. For M*=0.7 M Q and n oc i -2 , t > 10 3 yr, the mass 
of ionized gas increases with time. In these conditions, at 
each time step a new layer of molecular gas is exposed 
to the X-ray heating radiation, as H2 molecules are ad- 
vected from deep in the PDR slab toward the irradiated 
surface. Therefore, compared to the predictions of equi- 
librium calculations, a larger amount of hot molecular gas 
is formed, which emits stronger H2 vibrationally excited 
lines. This effect, discussed in detail in NH98, is larger in 



models with lower density, so that the 1-0S(1) intensity 
is higher in models with lower n. The opposite is true in 
models where the H 2 chemistry is treated under the as- 
sumption of stationary equilibrium. These models predict 
at t ~ 2500 yr a 1-0S(1) line about 7-10 times weaker (for 
no < 2.1 x 10 4 cm -3 ), mostly due to fluorescence in H2 
pumped by UV photons, and lower in models with lower 

The predicted intensities of all the H2 observed lines 
have been computed using a code which calculates the 
level population for the iV=299 bound states with rota- 
tional quantum number J <29 of the H2 molecule. The 
code include the effects of UV pumping by an external ra- 
diation field as well as collisions with H, H2, He, electrons 
and protons (Draine & Bertoldi 1996). We have used as in- 
put the physical conditions (namely, the radiation field at 
the inner edge of the PDR and the run with the depth in 
the PDR of temperature and fractional abundances of H, 
H2, He, electrons and protons) computed with the NH98 
code for t = 2500 yr. 

The models agree rather well with the observations 
for all the lines. Fig. |^ shows a Boltzmann plot for the 
H2 transitions, where the column density of the upper 
level of the transition (divided by its statistical weight) is 
plotted as function of the energy of the level. The filled 
dots are the observed values, the open symbols show the 
prediction of the best-fitting model (Mod 1, no = 7 x 10 3 
cm -3 ). The agreement is somewhat worse for models with 
higher density, which tend to predict weaker lines from the 
higher excitation levels. 



4.2. Shock Models 

In the NH98 models of the H2 shock emission, the relevant 
parameters that determine the intensity of the lines arc the 
shock velocity v s h (i.e., the difference between the torus 
expansion velocity and the red-giant wind velocity) and 
the pre-shock density, which, at any given distance R from 
the star, is determined by the red-giant wind properties 
(n w = M w /(4irR 2 iiv w f w ), where M w is the rate of mass- 
loss, v w the wind velocity, /,„ the fraction of solid angle 
over which the wind was ejected and /z the mean molecular 
weight). The H2 line intensity does not depend on the 
torus properties, but only on its expansion velocity, as long 
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Fig. 5. Boltzmann plot. The values of the column density 
of the upper level of the transition (divided by its statis- 
tical weight) are plotted as function of the energy of the 
level. The filled dots are the values derived from the ob- 
served lines. Arrows are 3cr upper limits. The prediction 
of the best-fitting model (no = 7 x 10 3 cm -3 ) are shown 
by the open triangles for lines of the v=0-0 band, open 
squares for the 1-0, open circles for the 2-1 and pentagons 
for the 3-2 band. 



as the matter in the torus does not become completely 
ionized or photodissociated. 

In NGC 2346, the expansion velocity of the neutral 
gas derived from CO data ranges from ~ 15 to >35 km 
s- 1 (Bachiller et al. 1989). We show in Fig. ||the predicted 
1-0S(1) intensity at t = 2500 yr as a function of the wind 
parameter M w / f w and two shock velocities, v s h=lQ km 
s- 1 (solid line), and w s /j = 15 km s~ 1 (dashed line). Higher 
values of v s h will appreciably dissociate H2 in J shocks. 
The intensity is computed according to Eq.17 of NH98. 
For v s h=lQ km s _1 , one needs M w /f w ~ 4.4 x 10~ 4 M 
yr _1 in order to reproduce the observed 1-0S(1) intensity. 
Unless f w <C 0.1, this implies an unusually high rate of 
mass- loss in the red-giant wind (Loup et al. 1993). The 
density of the pre-shock gas (n w ) is about 5 x 10 4 cm -3 ; 
this value is also unconfortably high, given the fact that 
the observed density of the CO torus (postshock gas), 
which these models predict to be ;§> n w , is <10 4 cm -3 
(Bachiller et al. 1989). A lower mass-loss rate is required 
if the shock velocity is higher. For v s h=15 km s — 1 , the ob- 
served intensity is reproduced by M w /f w ~ 1.7 x 10~ 4 M 
yr -1 and a pre-shock density of about 2 x 10 4 cm -3 . The 
shock models of NH98 predict a ratio of the 2-lS(l)/l- 
0S(1) intensity of 0.10 for v sh = 10 km b -1 , and 0.19 for 
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Fig. 6. Intensity of the 1-0S(1) line emitted in the shock 
between the expanding torus and the precursor red-giant 
wind at t = 2500 yr as function of the red-giant wind 
parameter M w /fw (see text). The top scale gives the cor- 
responding value of the pre-shock density n w for v w = 8 
km s _1 . The two curves correspond to two values of the 
shock velocity, as labelled. The observed 1-0 S(l) intensity 
is shown by the horizontal arrow. 



v s h = 15 km s , somewhat larger than the observed value 
(-0.07). 



5. Discussion 

The comparison of the observed intensity of the H2 T0S(1) 
line in the central torus of NGC 2346 to the NH98 model 
predictions shows that the emission can be produced in 
the hot PDR generated by the radiation of the central 
star, once the effect of X-ray heating and time-dependent 
(advecting) chemistry is taken into account. The best fit is 
obtained by models with relatively low density of the neu- 
tral gas (in agreement with the low density of the ionized 
material inferred by several authors). However, these PDR 
models require the luminosity of the hot central star to be 
significantly higher (L*— 250 Lq) than current estimates. 

If the H 2 lines are emitted in the PDR, we expect to 
observe a similar morphology in the ionized and H2 emit- 
ting gas. In NGC 2346, the PDR origin of the H 2 lines is 
supported by the fact that the same morphology is seen in 
H2 and in Ha (see Walsh 1983). Also, we detect Br7 emis- 
sion in the two H2 peaks, with a N-S profile that follows 
that of the 1-0S(1) line (Fig. |). The intensity of Br 7 pre- 
dicted by the models is very low, of the order of 8 x 10~ 6 
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erg cm 2 s 1 sr 1 , comparable to the observed values (4- 
6xl0~ 6 erg cm~ 2 s _1 sr -1 ). This supports our estimate 
ofL*. 

In principle, the observed intensity of the H2 1-0S(1) 
line can also be accounted for by the emission of the 
shocked gas produced by the expansion of the torus in- 
side a precursor red-giant wind. However, we estimate 
(following NH98) that one needs a rather high value of 
the mass-loss rate in the red-giant wind [M w j f w > 10 
Mq yr _1 ). This, in turn, implies a high density of the pre- 
shock gas (>10 4 cm -3 ), which is not supported by any 
existing observation. In fact, as discussed by Zuckerman 
and Gatley (1988), the main difficulty in ascribing the H 2 
vibrationally excited emission to shocks comes from the 
high momentum rate these models require. Many authors 
(see the review by Kwok 1993) have proposed that the for- 
mation and expansion of the PN shell (or torus) is related 
to the action of the fast wind from the central star. In this 
case, the ambient gas (pre-shock red giant wind) gains mo- 
mentum approximately at the rate at which momentum 
is delivered to the torus by the fast wind. Since the fast 
wind is radiation driven, this rate (P) must be <L*/c. In 
NGC 2346, we estimate that P is at least - 2 x 10 28 erg 
cm -1 s~ 2 , i.e., more than 600 times the present value of 
L*/c (for L*=250 L Q ) and 12 times higher than the max- 
imum L*/c reached by the star in its earlier evolution, 
according to the evolutionary tracks of Blocker (1995). 

The interpretation of the H2 emission in terms of shocks 
is often justified in the literature by the low measured ra- 
tio of the 2-lS(l) to the 1-0S(1) intensity. This argument, 
however, is not very strong, since in dense PDRs the low 
vibrational H2 levels are thermalized. The PDR models 
discussed in §4.1 predict a ratio 2-lS(l)/l-0S(l)~0.15 (not 
very different from the 0.10-0.18 range predicted by shock 
models; see §4.2), independently of the density and stel- 
lar luminosity. These values are somewhat higher than the 
observed ratio (~0.07). It is possible, and worth further 
investigations, that models tend to overestimate the fluo- 
rescent component of high vibrational lines, possibly be- 
cause of uncertainties in the collisional deexcitation rates. 

An interesting result of our observations, and one that 
we cannot account for with our simple models, is the vari- 
ation of the 2-lS(l)/l-0S(l) ratio with position along the 
slits, ranging from about 0.08 to 0.15 along the W slit, 
and from 0.08 to 0.23 along the E slit (see Fig. |). These 
variations are not monotonic with the distance from the 
peaks, but show evidence of structures, especially along 
the W slit. It is possible that this is due to density vari- 
ations, which affect the fluorescent contribution to the 2- 
1S(1) line. If the emission is due to shocks, this could trace 
variations in the propagation velocity of the shock in an 
inhomogeneous medium. 

Further support to the PDR origin of the H2 emission 
can be obtained by observing lines from higher vibrational 
states. We show in Table 3 PDR model-predicted values 
for lines not detected so far in NGC 2346, which, however, 



Table 3. High-v H 2 PDR Predicted Line Ratios 



Line A Mod 1 Mod 2 Mod 3 Mod5 



H (i-os(i)=ioo) 



4- 3S(l) 

5- 4S(l) 

6- 5S(l) 



2.5414 


15 


20 


20 


18 


2.7172 


0.6 


0.8 


0.8 


0.7 


2.9207 


0.2 


0.3 


0.3 


0.3 



Table 4. Mid-Infrared H 2 PDR Predicted Line Intensities 



Line 


A 


Mod 1 


Mod 2 


Mod 3 


Mod5 




H 


(10- 


5 erg cm" 


2 s" 1 sr" 




0-0S(0) 


28.22 


0.07 


0.06 


0.03 


0.05 


0-0S(l) 


17.03 


1.4 


1.3 


0.6 


1.2 


0-0S(2) 


12.28 


0.8 


1.0 


0.6 


0.7 


0-0S(3) 


9.66 


2.8 


3.1 


2.5 


1.9 


0-0S(4) 


8.02 


1.3 


1.0 


0.9 


0.5 


0-0S(5) 


6.90 


5.0 


2.7 


1.9 


1.0 


0-0S(6) 


6.11 


1.7 


0.7 


0.4 


0.3 



are accessible from space. The lines are very weak with 
respect to the 1-0S(1), with ratios that do not depend sig- 
nificantly on the density. We expect that they will be sub- 
stantially weaker in shocks. To complete the discussion of 
the H 2 spectrum, we show in Table 4 model-predicted val- 
ues of the intensity of mid-infrared lines in the v=0-0 band 
for four of the PDR models described in §4.1. These lines 
have been observed by ISO in a number of PNe (among 
them NGC 2346; Barlow et al. in preparation), and may 
be useful diagnostic of the physical conditions (see NH98). 

6. Summary and Conclusions 

This paper is motivated by our interest in understand- 
ing the origin of the H2 emission in PNe. Recently, Natta 
& Hollenbach (1998) have computed the evolution with 
time of the H 2 emission expected in PNe of different core 
mass and shell (or torus) properties. Their models com- 
pute the emission originating in the photodissociation re- 
gion (PDR) created by the UV radiation of the central core 
incident on the inner edge of the neutral shell, as well as 
the emission in the shock associated to the expansion of 
the torus inside the precursor red-giant wind. NH98 show 
that both regions can produce intense H2 emission and 
that detailed studies of individual PNe are necessary. 

NGC 2346 is a good candidate. It is a bright, young 
PN, with a typical butterfly morphology characterised by 
well developed bipolar lobes of emission and a bright torus 
around the central star. The same morphology is seen in 
Ha, CO and in the H2 1-0S(1) line. There is a large amount 
of information in the literature, concerning the properties 
of the ionized region and the gas kinematic that can be 
used in our analysis. 
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We have collected near infrared observations of NGC 23 
The data include broad K band image, an image in the 
H 2 vibrationally excited 1-0S(1) line and slit spectra in 
the K band in three positions in the nebula. The images 
confirm the well-known NGC 2346 morphology, with a 
central, bright torus, surrounded by weaker emission with 
a typical butterfly shape. The K band spectra show 11 
H2 lines with excitation energies from 6150 to 12552 K. 
Profiles of the lines intensity along the slit show evidence 
of secondary condensations outside the torus midplanc. 

A comparison of our observations with the NH98 mod- 
els shows that PDR emission can account for the H2 obser- 
vations. This requires a low-density shell (n < 10 4 cm -3 ), 
in agreement with the low density measured in the ionized 
region. Note that steady-state models or models where the 
soft X-ray radiation from the star is ignored predict a 1- 
0S(1) intensity one order of magnitude lower than models 
where both these effects are included. PDR models of the 
H 2 emission needs a central star significantly more lumi- 
nous than estimated in the literature (250 L Q against 17- 
90 L Q ). However, we think that this is consistent with all 
the available data, including our own Br7 observations. 

It is unlikely that the H2 emission originates in the 
shock between the expanding shell and the precursor red 
giant wind. Shock models require a much larger momen- 
tum input to the torus than possible from the central star. 
We estimate the discrepancy to be at least a factor 600 for 
the current luminosity of the central star (taken to be 250 
L Q ) and a factor 12 if we consider the maximum luminos- 
ity ever reached by the star in its previous evolution. 

In conclusion, we have proved that the PDR origin for 
the H2 emission in NGC 2346 is likely and that the mod- 
els, even if very simple, can account for a number of the 
observed properties. We suggest some additional obser- 
vations of lines from higher excitation vibrational levels 
and in the 0-0 band, which may help in determining the 
physical conditions in the shell with higher accuracy. 
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